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A general analysis on Fermion pairing at finite temperature 
and density between different species with mismatched Fermi 
surfaces is presented. Very different from the temperature 
effect of BCS phase, the recently found breached pairing phase 
resulted from density difference of the two species lies in a 
region with calabash-like shape in the T — fJ, plane, and the 
most probable temperature for the new phase's creation is 
finite but not zero. 

PACS numbers: 74.20.-z, 03.75.Kk, 12.38.-t 

The Fermion pairing between different species with 
mismatched Fermi surfaces, which has already been dis- 
cussed many years ago in investigating superconducting 
metal or B state of liquid 3 He in a magnetic field [1,2], 
prompted recently new interest in both theoretical and 
experimental studies. In high energy physics the strong 
interacting matter at high baryon densities described by 
Quantum Chromodynamics (QCD) is a color supercon- 
ducting phase, in which quarks of different colors pair 
with each other due to attraction in color anti-triplet 
channel [3,4]. The corresponding physical systems may 
be found in compact stars [5]. In non-relativistic case 
the study is related to the challenging goal of observing 
the BCS transition in trapped Fermionic atoms [6-9], the 
electrons in solid distributed in two different bands [6,9], 
and the neutron-proton pairing in isospin asymmetric nu- 
clear matter [10]. 

When the two Fermions have the same Fermi sur- 
face, the pairing phenomena is the well-known BCS 
mechanism [11]: At low temperature but high density, 
Fermionic matter with a sharp and high Fermi surface 
is unstable under attractive interaction(even very weak), 
a BCS paired state is favored instead. However, under 
some physical constraints the two Fermions may come 
from different Fermi surfaces. In atom traps, there can 
be two different hyperfine states of the same atom( 6 Li or 
40 Jf) which serve as two attracting species of fermions, 
and even mixture of two different Fermionic atoms could 
be realized [6-9]. In QCD matter at moderate baryon 
density (below the critical Color-Flavor-Locking density) 
which is considered to be relevant for understanding neu- 
tron stars, the mass difference between light quarks (u 
and d) and strange quark can not be neglected, and there- 
fore they have different Fermi momenta [12]. When chi- 
ral symmetry restoration at finite density is considered 
[13], the mass difference becomes more significant, and 
the mismatch will be enhanced. The mismatch happens 
even for the pairing between the light quarks when charge 
neutrality is taken into account [14]. 



A well-known pairing mechanism of Fermions from 
mismatched Fermi surfaces is the LOFF state [1]. Differ- 
ent from a BCS pair which has zero total momentum, a 
LOFF pair has finite total momentum. An analogue of 
this idea in dense QCD matter is the crystalline color su- 
perconductivity [15]. Recently a spatially uniform mech- 
anism is proposed [6] which leads to a breached pairing 
superfluidity (BP) [9] with coexistence of pairing state 
at lower and higher Fermi surfaces and particle state in 
between, when the densities of the two species differ from 
each other. This new phase is argued to be more favored 
than the LOFF phase for certain region of parameters. 

It is well-known that with increasing temperature 
of the system the superfluidity gap in a BCS state 
drops down monotonously and the region of superflu- 
idity phase in the temperature-density plane is reduced 
monotonously too [16]. What is the temperature effect 
of a system in a BP state and what is the difference from 
that in a BCS state? In this letter, we present a gen- 
eral analysis on spatially uniform pairing state with mis- 
matched Fermi surfaces at finite temperature and den- 
sity. We will discuss the stability of BP state induced 
by mass difference of the two species, and calculate the 
temperature behavior of the superfluidity gap and the 
phase diagram for a system with two species of different 
densities. 

We start with a system containing two species of 
fermions represented by a and b. The interaction can 
be modelled by a four-fermion point coupling, which is 
appropriate for both trapped Fermionic atoms and dense 
quark systems [3,6-9]. Since our purpose is a general 
analysis for pairing phenomena, we neglect inner struc- 
tures like spin, isospin, flavor and color, which are im- 
portant and bring much abundance while are not central 
for pairing. We write down the following Hamiltonian, 
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where p and q are momenta of the species, a, b, and w 
are the annihilation and creation operators, the coupling 
constant g is positive to keep the interaction attractive, 
V is the system volume and in continuous limit we sim- 
ply replace y J2p~ with J y 2 ? 3 , and the effective particle 

c a,b 



energies e" 



are */p 2 + m 2 a b — /i a ,b in relativistic case and 



2 ^ - — fj, a ^ in non-relativistic case. The Fermi momenta 

Pp b determined by = are controlled by the particle 
masses m a ,b and the chemical potentials fi a ,b- Whether 



1 



Ha,b can be used as free parameters depends on the phys- 
ical system we discuss. When the densities of the species 
are restricted by some physical constraint like fixed over- 
all particle density or fixed relative particle density, fj, a 
are not fully free and should be adjusted to satisfy the 
constraint. In fact, it is the density constraint which 
brings nontrivial BP state, as pointed out in [6-9] and 
will be clearly shown below. 

We introduce in the light of Mean-Field an order pa- 
rameter A = y < atb^p > and make it real by 
proper choice of phase factors of the creation operators, 
which allows us to diagonalize the Hamiltonian into 
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after a Bogliubov transformation from Fermions a and 
b into quasi-Fermions A and B with annihilation and 
creation operators A,B,A\W and quasi energies 
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The dispersions for particles a, b and quasi-particles 
A, B in relativistic case are shown in Fig.(l). With the 
interest in mismatched Fermi momenta induced by mass 
difference and density difference, a remarkable feature is 
that Ep' B can cross zero between the particle Fermi mo- 
menta Pp '. A simple algebra calculation shows that the 
condition for Ep^ to cross zero is the constraint A < A c 
for the gap. In relativistic case with the approximation 
of pp b » m a fi, the critical gap is 
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in non-relativistic case one has exactly 
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For convenience, we introduce momentum interval I = 
{p\p~ < \p\ < p+} with p± the solutions of E A ^ = 0. 

Whether Ep [B) crosses zero along the momentum axis 
plays a crucial role in pairing with mismatched Fermi 
surfaces. Ep^ = means directly gapless excitation 
of quasi-particles in a superfluidity phase observed in 
[4,6-9,17]. We see from Fig.(l) large mismatch and small 
gap leave space for gapless phenomena. If there is no mis- 
match, there is no place for gapless excitation. 

In order to understand the coexistence of pairing state 
and particle state when the gapless excitation region I is 



not empty, we write down the ground state of the system 
in terms of the zero quasi-particle state \0p ^ > and one 
quasi-particle state |lp^ > at zero temperature, 




FIG. 1. Dispersions of particles (dashed lines) and 
quasi-particles (solid lines). The thick and thin solid lines 
correspond to A < A c and A > A c , respectively. The left 
panel is for Ep = 0, and the right panel is for Ep = 0. 
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where 9(x) is a step function. With the help of the Bogli- 
ubov transformation the ground state can be expressed 
in the space of the original Fermions a and 6, 



.p-sinBpaW^lO^Op > x 



\g >= Y[ [cosb 
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U{o(PF-p F )al+e(p F -p F )bl P )\o;,o P > . (8) 

The above state tells us clearly the breached pairing phe- 
nomena. If neither Ep nor Ep 3 crosses zero, namely with- 
out gapless excitation, the species a and b are symmetri- 
cally and partially occupied and paired, the system is in 
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a BCS pattern. If one can find roots p± from E p 
there exists a momentum interval I in which one specie is 
fully occupied and the other is fully empty. The pairing 
occurs only in the region around the lower Fermi surface 
and the region around the higher Fermi surface, the two 
regions are separated by the gapless excitation region /. 
This means a new phase — the breached pairing super- 
fluidity proposed in [9] . The observed gapless phenomena 
[4,6-8,17] all fall into the BP phase, which is character- 
ized by two universal features, pairing breached by single 
occupation and simultaneous gapless and gaped compo- 
nents. 

We now turn to thermodynamics. From the Hamilto- 
nian (2) in quasi-Fermion representation the grand po- 
tential is easily expressed as 
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With given as a function of temperature, chemical po- 
tentials and gap parameter, all other thermodynamical 
quantities can be immediately obtained and its minimum 
gives the gap equation for solving A 
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with the Fermi-Dirac distribution f(x) = eX /T +1 - 

To do numerical calculation and present quantitative 
results we consider in the following only relativistic case. 
The calculation can straightforwardly be extended to 
non-relativistic case. Since the four-Fermion interaction 
is not renormalizable, one must use a momentum cut- 
off A to avoid divergence. We choose A = 600 MeV in 
our numerical calculation. The mismatch of the Fermi 
surfaces can be induced by either mass difference or den- 
sity difference of the two species. We now solve the gap 
equation (10) at finite temperature in the two cases. 

Mass Difference — We consider a system with u, a ^ b = 
\x but m a 7^ nib to study the effect of mass difference. 
For given ra a ,mf,,/u and T one can define two couplings 
9ci < 9c2 , corresponding to the critical conditions to have 
a BP state and a BCS state, respectively, 
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FIG. 2. The potential Q as a function of A at given T, fi 
with m a = 5MeV,m b = 200MeV, and g = 50 GeV~ 2 . 

For g < g c i the only solution of (10) is A = which 
means a normal phase without pairing. For g > g C 2 there 
is a stable BCS solution and A = becomes unstable. 
When g is in the region g c \ < g < g C 2, there are three 
solutions of (10): a normal solution A = 0, a BCS so- 
lution A > A c and a BP solution A < A c . However, 
only the BCS solution is the global minimum of f2, the 
BP solution is unstable because it doesn't satisfy the sta- 
ble condition d 2 £l/dA 2 > 0, and the normal solution is 
metastable. This is shown clearly in Fig. (2). When the 
temperature increases, the normal solution becomes sta- 
ble and the BCS solution changes to be metastable. In 



any case the BP solution can not be stable. When T 
exceeds some critical value, the unstable BP state even 
disappears. Therefore, the stability analysis rules out the 
possibility to create a BP phase by mass difference. 

Density Difference — From the thermodynamic poten- 
tial (9) it is easy to obtain the particle densities 
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We consider a system with fixed relative density 
— — A . 



(13) 



To satisfy this density constraint, namely a charge con- 
servation such as electronic charge with a, b carrying 
charge number — A, +1, respectively, the chemical po- 
tentials should be adjusted to be /i a = /i — A5/i and 
/Ub = fi + 5[i 7 where fi corresponds to the total number 
density since n a + ni, = —dfl/d/j,. Unlike the case of mass 
difference where the chemical potentials are free and A 
is purely determined by the gap equation, here only u, 
is free and A and (5/z are coupled to each other through 
(10) and (13). 

For A = 1 the system can be in a BCS state, while for 
A ^ 1 there will be some particles left after the pairing, 
the system must be in a BP state if pairing occurs. We 
first calculate A as a function of u, for A = 2 at given tem- 
perature. We choose the parameters m a = m\, = 5 MeV, 
and g — 80 GeV~ 2 in the following numerical calcula- 
tion. Introducing mass difference in calculation brings 
only quantitative changes. The value of g needed to pro- 
duce a BP decreases fast with increasing momentum cut- 
off A, and taking into account all inner degrees of freedom 
like spin, isospin, color and flavor will reduce the coupling 
by order of magnitude. From Fig. (3) a BP phase exists in 
a relatively narrow region, compared with a BCS state. 
The explanation for the existence of the higher critical 
chemical potential is that when \i is high enough the mis- 
match of the two Fermi surfaces is too large for the two 
species to pair at given coupling. 

In a BCS state it is well-known that the gap and the 
region of superfluidity decrease monotonously with in- 
creasing temperature. However, the case in a BP state 
is very different. With increasing temperature, the gap 
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firstly goes up and then drops down, and correspondingly 
the region firstly expands and then contracts. The tem- 
perature related to the maximum gap and the maximum 
region is a finite value but not zero. We can understand 
this feature in the following way. When the two Fermi 
surfaces coincide, the temperature deforms and lowers 
the Fermi surface and then the gap decreases. In the case 
of mismatch, the two different sharp Fermi surfaces are 
deformed and lowered by the temperature on one hand, 
but the expansion of the two distributions make the two 
species much closer in the phase space and then much 
easier to pair on the other hand. The former temper- 
ature effect reduces the gap but the latter effect favors 
the gap. The competition between the former and the 
latter effects controls the temperature dependence of the 
amplitude and the region of superfluidity. 

The behavior of the gap as a function of temperature 
at fixed chemical potential for different relative density 
is shown in Fig. (4). At A = 1, the gap drops down 
monotonously due to disorder brought in by increasing 
temperature, a standard characteristic of BCS phase. For 
A ^ 1, we see again the amazing temperature behavior: 
The gap first increases to a maximum value and then re- 
duces rapidly to zero with increasing temperature. While 
the density difference A ^ 1 is the prerequisite for a BP 
phase, the amplitude and the region of the gap decreases 
with increasing A. 

We further give in Fig. (5) the phase diagram in T — \i 
plane for the BP phase with A = 2 and compare it with 
the BCS phase diagram with A = 1. The calabash- like 
phase transition line from normal state to BP state man- 
ifests the competition between the two opposite tempera- 
ture effects on the mismatched Fermi surfaces, discussed 
above. The most probable temperature for producing a 
BP state is not zero but a finite value. 

In summary, we have presented a general analysis 
on Fermion pairing between different species with mis- 
matched Fermi surfaces under the assumption of spatial 
uniform, and investigated especially the temperature be- 
havior of the new BP phase. From the stability analysis 
the mass difference can not create a stable BP state. In 
the case of fixed relative density rib/n a ^ 1 the tempera- 
ture effect not only deforms and reduces the mismatched 
Fermi surfaces which leads to the usual suppression of the 
gap, but also make the overlap region of the two species 
even wider which favors the condensate. The competi- 
tion of the two opposite temperature effects results in 
a calabash-like phase transition line from normal to BP 
states. The temperature corresponding to the largest gap 
and the largest region of superfluidity is finite but not 
zero. Due to this novel temperature effect, we think that 
trapping different Fermionic atoms with different densi- 
ties at controlled finite temperature may be a feasible 
way to obtain a BP phase. 
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FIG. 3. The gap parameter as a function of [x at fixed rel- 
ative density A = 2 for different temperature. 
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FIG. 4. The gap parameter as a function of T at fixed 
chemical potential for different relative density. 
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FIG. 5. The BCS (X = 1) and BP (X = 2) phase transition 
lines in T — (i plane. 
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